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Abstract 

We have studied the response of Tl-doped and Na-doped Csl crystals to nuclear 
recoils and 7's below 10 keV. The response of Csl crystals to nuclear recoil was 
studied with mono-energetic neutrons produced by the 3 H(p,n) 3 He reaction. This 
was compared to the response to Compton electrons scattered by 662 keV 7-ray. 
Pulse shape discrimination between the response to these 7's and nuclear recoils 
was studied, and quality factors were estimated. The quenching factors for nuclear 
recoils were derived for both CsI(Na) and CsI(Tl) crystals. 
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1 Introduction 



It is known that a major component of matter in the Universe is not ordinary 
matter(luminous stars and/or baryons) but nonbaryonic exotic matter[l]. The 
Weakly Interacting Massive Particle(WIMP) is one of the strongest candidates 
for this dark matter, and direct searches for WIMPs using various detection 
techniques have been undertaken [2-7]. 

The sensitivity of WIMP detection strongly depends on various parameters 
of the detector, such as detection threshold, mass of the detector, background 
rate, pulse shape discrimination(PSD) capabilities, etc. A WIMP may interact 
with the target nuclei through a WIMP-nucleus elastic scattering, and the 
detection of the nuclear recoil signal may give a signature of WIMP detection. 
The expected recoil energy is several tens of keV, depending on the masses 
of the target nucleus and the WIMP. However, because of the high stopping 
power of the recoil nucleus, the light yield of a scintillation detector is 5-10 
times lower in comparison to a signal from a 7 of equivalent energy due to 
quenching effect in normal crystals. Thus, the detection threshold should be 
low enough to measure a few keV of energy. 

The expected event rate of the WIMP-nucleus interaction is less than 1 event 
/kg/day. Such a low event rate requires an extremely low background rate, 
on the order of several counts/kg/day/keV, and a large detector size in order 
to increase the signal rate. To reduce the background rate, the PSD between 
nuclear recoils and background 7's or electrons in the several keV region, can 
be utilized. 

A Csl crystal is a good candidate for a WIMP detector thanks to its large 
light yield and relatively good PSD capability. Also, it is relatively easy to 
fabricate large detector volumes because the crystals are not hygroscopic. The 
large masses of Cs and I nuclei enhance the cross-section for a WIMP-nucleus, 
spin-independent interaction (which depends on A 2 [8]). However, the detailed 
response of the crystal to low energy radiation has not been fully investigated. 
Recently, the feasibility of using CsI(Tl) crystals as a WIMP detector was 
shown[9], while several other groups reported the characteristics of a CsI(Tl) 
crystal for WIMP detection[10,ll]. 

We studied here the response of Csl crystals with several different doping 
concentrations of Tl and Na to neutrons and 7's. The energy region was several 
keV for the 7's and several tens keV for the nuclear recoils. In this paper, the 
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measurements of PSD and the quenching effect of the nuclear recoils for both 
CsI(Tl) and CsI(Na) crystals are reported. 



2 Experiments 

The response of Csl crystals to neutrons was studied using a tandem accelera- 
tor(maximum terminal voltage of 1.7 MV) at the Korea Institute of Geoscience 
and Mineral resources(KIGAM). Mono-energetic neutrons of 2.62 MeV were 
produced at zero degrees via the 3 H(p,n) 3 He reaction, using a 7 to 9 nA beam 
of 3.4 MeV protons on a 3 H target. These neutrons produce nuclear recoils 
with energies near that expected for WIMP interactions. 

Figure 1 shows a schematic of the experimental setup. A Csl crystal was 
located at zero degrees with respect to the proton beam line. The Csl crystal 
was 3 cm x 3 cm x 3 cm, in cubic shape, and with all 6 faces polished. 
Two 3- inch PMTs(D726UK, Electron Tubes Limited), with green-extended 
RbCs bi-alkali photocathodes, were directly attached to the top and bottom 
surfaces. The top and bottom sides of the crystal were completely covered 
by the PMT cathode plane, and the other four sides were wrapped with two 
layers of Teflon (0.2 mm thick) followed by black vinyl sheets. In order to 
identify neutrons scattered from Csl, six neutron detectors were located at 
various angles ranging from 45 to 130 degrees. The detailed locations of the 
neutron detectors are shown in Figure 1. The neutron detectors were made of 
liquid scintillator(BC50lA) contained in a cylindrical Teflon container of 12.7 
cm diameter x 10 cm long. One 2-inch PMT(H1161, Hamamatsu Photonics) 
was used for readout. Each neutron detector was contained in a polyethylene 
block, about 5 cm thick, and was shielded with lead and borated paraffin 
blocks. The recoil energy (E recoi i) of the nucleus can be calculated by a simple 
kinematical equation using the incident neutron energy and the scattering 
angle of the neutron, as determined from the position of neutron detector: 

m n + m N J j ' 

where E beam is neutron energy, m n and m N are the masses of neutron and 
recoiling nucleus(Cs or I)Q respectively, and 9 is the neutron scattering angle. 

The signals from a Csl crystal are amplified using home-made, fast amplifiers 
(x8) with low noise and high slew rates. One output from these amplifiers 
(x8) was directly connected to a digital oscilloscope(LeCroy LT364) in order 

1 In the real calculation, the average mass of Cs and I was used. 
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to record the pulse shape. The other signal is amplified (xlO), again using a 
Philips amplifier(Philips 770), and connected to a discriminator(Philips 711) 
to form the trigger logic. The discriminator threshold was set to the level of 
single photoelectrons. In order to suppress the accidental trigger from PMT 
dark current, at least 4 photoelectrons within 1 /xs( actually 2 photoelectrons 
within 0.5 /is for each PMT) were required(a so-called "4-fold" coincidence). 
The detailed description is shown in the reference[9]. 

The event rate of 4-fold coincidences of the signals from the crystal was about 
1200 Hz when the neutron beam was on. Most of these events were due to the 7 
backgrounds induced by the neutron beam. In order to reduce the trigger rate, 
a time coincidence of a Csl signal (4-fold) and the neutron detector signal must 
occur within 500 ns. This reduced the trigger rate to less than 3 Hz. Trigger 
timing was determined by the neutron detector. The Csl signal was digitized 
by the digital oscilloscope(LeCroy LT364) with a 100 MHz sampling frequency 
and a 10 /is full range. The trigger timing for the oscilloscope was set at 3 
/is. The neutron detector signal was digitized by a charge sensitive camac 
ADC(LeCroy 2249A). Data acquisition used ROOT[13] on a linux machine. 

In order to compare the pulse shape of a nuclear recoil to that of a 7, the 
crystal was irradiated with a 137 Cs source to induce Compton scattering. This 
induced a signal which is distributed uniformly inside the crystal. These data 
were taken separately using the same setup, but without the proton beam. 
The trigger was a 4-fold coincidence from a crystal. 

Table 1 shows a list of the crystals tested. Two doping materials, Na and Tl, 
with various concentrations, were tested. The crystals were supplied by the 
Institute for Single Crystals in Kharkov, Ukraine. Growing methods of each 
crystal are indicated in the table. 

The doping concentrations of the CsI(Tl) crystals were measured by the ICP- 
AES (Inductively Coupled Plasma - Atomic Emission Spectrometry) method. 
The 4 corners of each crystal were scratched by a razor blade, and sent to two 
independent institutes. Each institute measured two samples for each crystal. 
The errors of the doping concentration shown in the table are the standard 
deviation of the 4 measurements. 

For CsI(Na), the doping concentrations were measured by ICP-MASS (Induc- 
tively Coupled Plasma - Mass Spectrometry) method. The measurement was 
done by only one institute and the measurement errors were claimed to be 2-3 
%. 



3 Analysis 
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3. 1 Typical low energy signal from a Csl crystal 

The single photoelectron(SPE) signal from a PMT has less than a 10 ns full 
width. After amplification, the signal is stretched a bit due to the inherent 
shaping time of the amplifier(the rise time is about 10 ns and the fall time 
is about 20 ns). Hence, the SPE signal spreads about 30-40 ns in the time 
spectrum. It is known that the decay time of CsI(Tl) is about 1 /xs [14]. Since 
the number of photoelectrons for low energy events (less than 10 keV where 
WIMP signals are expected to occur) was less than 60 (or 30 per one PMT), 
the SPEs were clearly separated from each other in an event. That is, a low 
energy real event can be identified as a set of SPEs above a threshold cut 
within a timing window of 10 /is. 

An example of a low energy signal is shown in Figure 2. This signal was 
produced by 5.9 keV X-ray from a 55 Fe source. The average number of photo- 
electrons from a CsI(Tl) crystal was about 5 p.e./keV. As expected, the signal 
was composed of a set of well separated peaks composing the SPE signal. In 
order to reduce the noise contribution, a 'clustering algorithm' was developed 
to identify SPEs. The pedestal of each event was estimated in the first 2 /j,s 
window before the trigger for each PMT. Each peak of about a 30-40 ns width 
and above a threshold value was considered one cluster. 

Figure 3 shows the pulse height distribution of the clusters composed of hav- 
ing 4 keV < E meas < 8 keV, where E meas is the measured electron-equivalent 
energy. The figure shows separation between the SPEs and the noise, with a 
peak to valley ratio of about 2. If the cluster is a real SPE, the pulse height dis- 
tribution should follow a Poisson distribution, because the amplification at the 
first dynode inside the PMT follows Poisson statistics. In an attempt to under- 
stand the shape of the spectra, a fit was done using two superimposed Poisson 
functions(one for the SPE and the other for the SPE-overlapped signal) and 
an exponential function for noise, as expressed in the following equation. 

f = A -w^ + B -w^s + c - e ' x '^ 

T(r + 1) T(r + 1) 

r = x ■ g,n = m- g,fi' = rri ■ g, 

where x is the pulse height, m{m!) is a number corresponding to the mean 
of the Poisson distribution, and g is the gain factor of the PMT (from the 
second dynode). The results of the fit are overlaid in the figure as a solid line, 
where the reduced x 2 °f the fit is about 1. The ratio of the mean values of 
two Poisson distributions, m/m', was 1.88, compared to an expected value 
of 2. The ratio of the contributions of the two Poisson distributions is about 
12%. So, we could conclude that the cluster is a SPE (~ 80%) or due to 
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overlapped photoelectrons (~ 20%). The same fit was performed for events of 
2 keV < E meas < 4 keV. A single Poisson function will fit the SPE spectrum. 

3.2 Neutron identification 

In order to identify neutron induced events, we required a coincidence between 
one of the 6 neutron detectors and the Csl crystal. Furthermore, neutron scat- 
tering events can be selected by separating neutrons and 7's in the neutron 
detector. It is well known that liquid scintillator(BC501A) can separate neu- 
trons from 7's using PSD. The neutron (proton recoil) events have a longer 
tail than 7 events for the same energy loss[14] in the BC501A liquid scintil- 
lator. Therefore, using the ratio of the total charge to the charge in the tail 
of the signal, the neutrons can be clearly separated from 7's. Figure 4 shows 
the neutron identification plot for the neutron detector. Above 300 keV, the 
neutron-7 separation is about 2.7a which is sufficient for neutron-event se- 
lection. The detailed description for neutron detector will be described in a 
separate report [15]. 

Figure 5 shows the recoil energy spectra of the CsI(Tl-0.073mole%) crystal 
tagged by a neutron detector. The measured energy increases as the scattering 
angle increases as expected. A quantitative analysis will be described in Section 
5. 

3. 3 Timing characteristics of Csl crystal 

To study the overall timing characteristics of a crystal at low energy, we accu- 
mulated whole events for each crystal. Figure 6 shows the accumulated time 
spectrum of a CsI(Tl-0.073mole%) crystal for neutron induced events, whose 
measured energy is 4 keV < E meas < 10 keV. In the figure, the time scale is 
in fis, and the timing of the neutron detector was set to 3 jis. 

The timing characteristics of the Csl response can be parameterized by a 
rise time(r r ), two decay times(r/ and r s ), and the ratio of slow to fast decay 
components(-R): 

A f 1 R ^ f 1 R \ *-*0 I , 

F(t)=A-l — -e T f + — -e —-[ — + — ] -e —\ + b.g., (3) 

where time zero(to), overall normalization(A), and the amount of background 
photoelectrons (fe.g.) are included. The amount of background photoelectrons 
was estimated from the number of photoelectrons within a 2 /is window during 
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the 1 /is prior to the event. The constant b.g. values were obtained by averaging 
all the events under analysis, resulting in about 0.1-0.3 photoelectrons per 1 
/is. The fit was over 10 fis even though the figure shows only 5 /xs. The values 
of parameters from the fit are shown in Table 2. For the Tl-doped crystals, the 
rise time is about 20 ns and it has two components of decay times, fast (about 
0.54 /is) and slow(about 2 /is) components. For the Na-doped crystals, the rise 
time is about 30 ns, the fast decay component is about 0.41 /is, and the slow 
decay component about 2.4 /is. The time characteristics of CsI(Na) are similar 
to those of CsI(Tl) at low energy. This is not an expected feature because it 
is well known that the decay time of CsI(Na) is much shorter than that of 
CsI(Tl) at high energy [14]. 

We do not observe a significant dependence on the doping concentration except 
for the case of CsI(Tl-0.071mole%). Both CsI(Tl-0.071mole%) and CsI(Tl- 
0.073mole%) have the same doping material and the same concentration, 
and the only known difference is the growing method, Kyropoulos(CsI(Tl- 
0.071mole%)) and Bridgman(CsI(Tl-0.073mole%)). However, it is not clear if 
different doping methods cause different time characteristics because we tested 
only one CsI(Tl) crystal made by Kyropoulos method. 



4 Pulse Shape Discrimination 

4-1 Comparison between nuclear recoil and 7 signals 

Pulse shape discrimination(PSD) with a Csl crystal is possible due to dif- 
ferent timing characteristics between neutron induced (nuclear recoil) and 7 
induced events[16]. Figure 7 shows the accumulated time spectrum of a CsI(Tl- 
0.073mole%) crystal for 4 keV < E meas < 10 keV, where the time spectrum for 
7 events and that for nuclear recoil events are overlaid. For each event, the tim- 
ing of the first photoelectron was used as to- The ambiguity due to improper 
determination of to should affect the time distribution, especially for the low 
energy event which does not have a sufficient number of photoelectrons. Eq. 
4 shows the fit function for the time spectrum, which is similar with Eq. 3. In 
order to reduce the bias on the t determination, the fit was done over 0.3 /is. 
Each component of the time spectrum is drawn as a solid line(nuclear recoil) 
or a dashed line(7). The amount of background photoelectrons for both cases 
are also shown. The fit shows that nuclear recoil events clearly have a shorter 
tail than 7 events. 

/ 1 -±- R t\ 
F(t) = A ■ — ■ e T f + - • e~ + b.g. (4) 

\ r / r s J 



7 



4-2 Quality factor 



To quantify the capability for neutron-7 discrimination among the different 
crystals, and to evaluate the sensitivity of a WIMP search, a quality factor 
was defined [17] as 



where a is the fraction of signal events passing the event selection criteria and 
f3 is the fraction of background events which passed the same criteria. For an 
ideal detector, a = 1 and (3 = 0. Therefore, a smaller quality factor means a 
better separation between signal and background events. 

For a simple separation parameter, a mean time is defined as; 

<*>=^-*o» ( 6 ) 



where Ai is the charge of the i-th cluster, U is the time of the i-th cluster, and to 
is the time of the first cluster (assumed as time zero). Figure 8 shows the mean 
time distribution for nuclear recoil and 7 events for a CsI(Tl-0.073mole%) 
with events between 6 keV < E meas < 7 keV. Nuclear recoil and 7 events 
show different mean times. By selecting the proper signal window, we can 
choose the event sample which maximizes the signal to background ratio. 

However, due to background photoelectrons which are uniformly distributed, 
the mean value of the < t > distribution is longer than the real mean time. 
To estimate the effect of background photoelectrons, the mean value of the 
mean-time distribution is calculated with the time constants(r/, r s , and R) 
obtained from the fit shown in Eq. 4 for each energy bin. The variance of 
the mean time distribution is calculated by removing the contribution from 
the background photoelectrons by assuming a quadratic sum(cr| t> = cr^ ignal + 
o~l g ), where cr% t> , o~b. g ., an d °~ signal are the variance of the < t > distribution, 
the background photoelectrons, and the signals, respectively. Figure 9 shows 
the mean time of the nuclear recoil and 7 events calculated from the time 
constant of the crystal. The errors are only statistical. 

Figure 10 shows the quality factors for CsI(Tl) for various doping concen- 
trations. The results from previous measurements for CsI(Tl)[10] and for 
NaI(Tl)[18] are overlaid. For CsI(Tl), the present results show a bit more 
energy dependence than Ref. [10]. However, in the energy region less than 5 
keV where most of the WIMP signals are expected, both measurements are 
comparable. CsI(Tl) shows in comparison to Nal(Tl) about a 10 times smaller 
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quality factor. In our measurement range, we do not observe a significant 
dependence of the quality factor on doping concentrations. 

Figure 11 shows quality factors for a CsI(Na) crystal. The quality factor of 
CsI(Na) is about 2 times larger than that of CsI(Tl), but still about 5-10 times 
smaller than that of Nal(Tl). The quality factor for a CsI(Na) crystal is about 
2-3 times worse than CsI(Tl) crystal in the E meas < ~7 keV region, and is 
almost comparable above that region. 



5 Quenching factor 



The light yield as a function of recoil energy (E recoi i) for various crystals is 
shown in Figure 12. The recoil energy is calculated by Eq. 1. As expected, 
CsI(Na) shows about a 30% larger light yield than CsI(Tl) for the same recoil 
energy. However, the light yield of the crystals with different doping con- 
centrations are consistent with each other within errors, except for CsI(Tl- 
0.426mole%), whose light transmission was not good. It seems that the light 
yield is already saturated as indicated by the doping concentration in our 
measurement range. 

The energy calibration to obtain E meas uses a 5.9 keV X-ray from 55 Fe for 
the CsI(Tl) crystal. Non-linearity of the light yield of the CsI(Tl) crystal at 
low energy[19] is considered as a systematic error of E meas . For the CsI(Na) 
crystal, the energy calibration uses a 59.5 keV X-ray from 241 Am since CsI(Na) 
is not sensitive to low energy photons, due perhaps to some unknown surface 
effect. The response of CsI(Na) is assumed to be linear within our measured 
energy range because the low energy response(_E meas < 10 keV), including the 
non-linearity for CsI(Na), has not yet been studied. 

Using E meas and E recoih the quenching factor is calculated to be; 



Figure 13 shows that our measured quenching factors for CsI(Tl) are in good 
agreement with previous measurements[10,12] within errors. The quenching 
factors for a Na-doped Csl crystal, which were measured for the first time, 
are shown in Figure 13. For the comparison, the quenching factors for CsI(Tl- 
0.071mole%) are shown together. The quenching factor of CsI(Na) is about 
50% smaller than that of CsI(Tl). 
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6 Conclusion 



The intensive studies on the response of CsI(Tl) and CsI(Na) crystals were 
performed at low energies (E meas < 10 keV), and the low energy response for 
a CsI(Na) crystal is reported for the first time. The pulse shape discrimina- 
tion(PSD) between 7 and nuclear recoil events was measured and the quality 
factor was estimated using PSD. The PSD and quality factor for CsI(Tl) are 
similar to previous measurements. The quality factor for a CsI(Na) crystal is 
about 2-3 times worse than for a CsI(Tl) crystal in E meas < ~7 keV region, 
and almost comparable above that region. 

The quenching factor for CsI(Tl) and CsI(Na) crystals of various doping con- 
centrations were also measured. As expected, CsI(Na) shows larger light out- 
put than CsI(Tl) even in the low energy region. However, the quenching fac- 
tor of CsI(Na) for nuclear recoil is about a factor of two smaller than that of 
CsI(Tl). 
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Table 1 



List of Csl crystals tested 


Doping Material 


Concentration (mole% ) 


Growing method 


Tl 


0.073±0.02 


Bridgman 


Tl 


0.128±0.01 


Bridgman 


Tl 


0.426±0.23 


Bridgman 


Tl 


0.071±0.01 


Kyropoulos 


Na 


0.0188 


Kyropoulos 


Na 


0.0254 


Kyropoulos 


Na 


0.0258 


Kyropoulos 
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Table 2 

Timing characteristics of Csl crystals for nuclear recoil(Errors are only statistical). 



Doping(mole%) T r (fj,s) T f(^ s ) t s ([xs) R mean time(^s) 



Tl(0.073±0.02) 
Tl(0.128±0.01) 
Tl(0.426±0.23) 
Tl(0.071±0.01) 
Na(0.0188) 
Na(0.0254) 
Na(0.0258) 



0.021 ±0.001 
0.017 ±0.001 
0.014 ± 0.002 
0.017 ±0.001 
0.035 ± 0.002 
0.029 ± 0.001 
0.032 ± 0.001 



0.546 ± 0.007 
0.560 ± 0.008 
0.523 ±0.020 
0.457 ±0.017 
0.432 ± 0.006 
0.409 ± 0.005 
0.433 ± 0.003 



2.193 ± 0.054 
2.240 ± 0.067 
1.653 ±0.082 
1.341 ±0.048 
2.437 ±0.050 
2.422 ± 0.046 
2.753 ± 0.036 



0.426 ±0.011 
0.393 ±0.012 
0.582 ± 0.024 
0.763 ± 0.021 
0.521 ±0.008 
0.463 ± 0.007 
0.474 ± 0.005 



0.981 ±0.011 
0.969 ±0.012 
0.924 ±0.025 
0.849 ± 0.020 

1.037 ±0.018 
0.966 ±0.012 

1.038 ±0.013 



13 



40keV 




Fig. 1. Experimental Setup. 
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from 55 Fe source. 
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Fig. 3. Single photoelectron spectrum: The fit results are overlaid(solid curve), done 
using two Poisson functions. 
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Fig. 6. Time spectrum from CsI(Tl-0.073mole%). The relative time of each event is 
fixed by one neutron detector at 3 fis. The solid curve represents the fit result by 
Eq. 3. 
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Fig. 7. Time spectra from CsI(Tl-0.073mole%) for nuclear recoil and 7 events. The 
timing of the first photoelectron of each event is fixed at 0. Each decay component 
and the background are decomposed by a fit to Eq. 4. The solid lines are for nuclear 
recoil data and the dashed lines are for 7 data, where the background for 7 events 
is less than 0.1. For a better display, only one out of every 5 points in the X-axis is 
shown in the figure. 
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Fig. 8. Mean time distribution for CsI(Tl-0.073mole%) for nuclear recoil and 7 
events. The mean time distributions are fitted well with the Log-Gaussian func- 
tion(solid curves). 
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Fig. 9. Mean time of CsI(Tl-0.071mole%) for each energy bin, calculated from the 
time constant of each crystal. The errors in Y-axis are statistical only. The error 
bars in X-axis show the energy region used for each data point. 
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Fig. 10. Quality factors for various CsI(Tl). The errors are only statistical. The 
present results(open markers) are compared to the data of S. Pecourt et al. [10] and 
of G. Gerbier et al. [18]. 
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Fig. 11. Quality factors for various CsI(Na). The errors are only statistical. 
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Fig. 12. Light yield vs recoil energy for CsI(Tl) and CsI(Na) with different doping 
concentrations indicated in the figure. 
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Fig. 13. Quenching factors for CsI(Tl) (present results: closed markers). Previous 
measurements (open markers) are compared. 
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Fig. 14. Quenching factors for CsI(Na) 
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